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1. INTRODUCTION 

Livestock farming is an activity of keeping animals that can be utilized by humans mainly as a food 
producer [1]. Two among many of the farm product is meat and eggs. In Indonesia, such consumption of 
meat products has risen consistently over the years. Among the meat product, 62% of meat product coming 
from poultry farm [2]. This statement supported by data published by USAID whose state that in 2013, 
Indonesia’s poultry industry produced about 1.285 million tons meat for a year. Meanwhile the production of 
poultry eggs in 2013 reached 1.2 million tons [3]. 

When raising livestock such as chicken or cow, one aspect that must be considered is the aspect of 
feeding livestock. The feed must have sufficient nutrition to keep livestock grow well, but it requires 
considerable daily cost, which is about 65-70% of operating costs [1], [4]. Based on the fact, farmers have to 
figure out how to meet the nutritional needs of cattle at the minimum cost. As farmer and expert know, 
animal feed usually are a mixture of feed materials from agricultural and industrial waste, if necessary 
additional components such as vitamins and mineral extracts may be added. The main nutrient content that 
must be considered is the protein, fat, and fiber [5]. But to make a good feed composition we also need to 
consider about some factors such as fluctuating prices, and proper nutrition simultaneously [4]. 

Optimal feeding composition is being needed to gain weight for both broiler chicken and laying 
hens. For broiler chicken, the more their weight gained the more carcass’ meat they produced, meanwhile for 
laying hens the more they gained weight, they will produce bigger egg and with better quality. For laying 
hens, after they reached optimal weight, the feed will act as energy provider for egg production.[6], [7]. 
While it possible to use manual search to mix feed ingredients, study already proved that the manual search 
of a solution for optimizing requires a great deal of insight and patience [8]. To solve this problem, some 
methods are already be introduced. The most classical approach to solve feed optimization problem usually 
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using Pearson’s Square (PS). PS is relatively easy to understand and implemented but it usually only 
maximizing the requirements for one or two elements of nutrition only. Besides if we judge from the price 
point, the price of the material sometimes becomes quite expensive [9]. 

Another proposed solution is using Evolution Strategies (ES). ES is part of evolutionary algorithms, 
a group of algorithm that has been successful to solve many composition problems on many disciplines such 
as Engineering, Biomedic, Economy, Operation Research, Social Science, Physics and many more. For 
example for drugs composition, livestock feed composition and Cutting Stock Problem [10]-[12]. ES itself 
was being invented by Ingo Rechenberg from Technical University Berlin at almost the same time with 
Genetic Algorithm’s (GA) invention [13]. While GA is widely used by many researchsrs, study by Salimans 
et.al. mention that ES may be used to be alternative evolutionary algorithm besides GA [14]. 

ES algorithm aims to find the optimum solution and can be applied as a solution to solve complex 
problems with many parameters and restrictions [1], [15]. Some research about this algorithm to solve 
combinatorial problem was already done such as work by Ahire et al [16]. which using ES to solve 
Workforce-constrained Preventive Maintenance Scheduling that usually exist on heavy equipment overhaul 
facilities such as aircraft service centers or railroad yards. On the other side, Schmidt-Dannert et al. [17] was 
trying to implement ES method to solve a very complex problem : create carotenoid diversity in E. coli cells. 

Another nature inspired algorithm besides ES that can be used is Firefly algorithm (FA). FA was 
relatively new among nature inspired algorithm. It was being invented by Xin-She Yang from Cambridge 
University on 2007. By took concept of fireflies flashing light that will become weaker after some distances, 
the fireflies object that represents real world problems will be attract into nearest firefly object with strongest 
flashing light [18]. FA also has many similarities with particle swarm optimization (PSO) algorithm, because 
it is also based on swarm behavior but FA has better efficiency and robustness in complex optimization 
problems [19]. On previous study by Wang et.al, FA with some modified attraction rule proved to be able to 
solve optimization function such as Sphere, Rosenbrock, and Griewank with relatively high success rate [20]. 

Study by Genfu et.al proved that FA can be used to solve NP hard problem, especially—in this 
study—Reactive Power Optimization. Meanwhile, research by Apostolopoulos and Vlachos [21] that trying to 
solve energy problems by solve the economic emissions load dispatch problem using the FA proved that FA 
is able to generate optimal solution with proper selection of control parameters. This study has objective to 
compare performance of ES and FA based on their average fitness and computational time. The average 
fitness test will measure how good the solution resulted by the algorithm while computational time will 
measure how fast the algorithm will give the result. 


2. RESEARCH METHOD 
2.1. Nutrient requirements 

To ensure the nutrient requirements for broiler and layer chickens can be fulfilled, we use some 
standard that can be viewed on Table 1 and Table 2. for broiler chicken case. For laying hens that produce 
eggs, the standard can be viewed on Table 3. [22], [23]. 

To make chicken feed we will combine four ingredients as can be seen on Table 4. It consists of 
corn, sorgum, rice, and coconut oilcake. The price is on Indonesian Rupiahs (IDR) and the price will be 
counted as price per kilograms. 


Table 1. Broiler Chicken Nutrients Requirement Table 2. Broiler Chicken Nutrients Requirement 
(Starter) (Finisher) 

No. Parameter Measure Constraint No. Parameter Measure Constraint 

1 Water % Max 14.0 1 Water % Max 14.0 

2 Raw Protein % Min 19.0 2 Raw Protein % Min 18.0 

3 Raw Fat % Max 7.4 3 Raw Fat % Max 8.0 

4 Raw Fiber % Max 6.0 4 Raw Fiber % Max 6.0 

5 Ash / Minerals % Max 8.0 5 Ash / Minerals % Max 8.0 

6 Calcium % 0.90 — 1.20 6 Calcium % 0.90 — 1.20 

T Total Phosphor % 0.60 — 1.00 7 Total Phosphor % 0.60 — 1.00 

8 Phosphor % Min 0.40 8 Phosphor % Min 0.40 

9 Metabolism Energy kkal/kg Min 2,900 10 Metabolism Energy kkal/kg Min 2,900 
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Table 3. Layer Chicken Nutrient Requirements 


Nutrient Requirements Phase 

Starter Layer Layer Finisher 
Water (%) <14 <14 <14 
Metabolism Energy (kkal/kg) 2,980 2,750 2,750 
Raw Protein (%) 20.5 18 17 
Raw Fiber (%) <6.5 <7 <7 
Raw Fat (%) 25-7 3.5-3.7 2.5-7 
Ash (%) 5-8 10- 14 10-14 
Calcium (%) 1.05-1.1 3.25-4 3.6-4 
Phospor (%) 0.48 0.40 0.33-0.37 


Table 4. The Ingredients 
Ingredients Protein Fat Fiber Ca P EM Price (IDR) 


Corn 8,6 3,9 2 0,09 0,37 3370 3000 
Sorgum 10 2,8 2 0,41 2,13 3250 2500 
Rice 8 1,7 9 0,09 0,08 2660 3000 
Coconut oilcake 14,1 8,12 23,32 0,29 0,79 1450 4000 


2.2. Representation 

The first step in implementing evolutionary algorithm is to define the exact representation of 
chromosomes (encoding). A good chromosomes representation is very important because it will affect the 
effectivity of ES and FA in exploring the search space [24]. To formulate the feed composition, we will use 
real code representation that describes the amount of each feed ingredient used in the feed mixing. 
The number of genes in a chromosome both in ES or FA is the percentage of feed ingredients and each 
chromosome will represent how many that ingredients will be used. Zero means that ingredients will not be 
included. The example of chromosomes can be seen on Table 5. 


Table 5. Gen Representation 
Corn Sorgum Rice Oilcake 
0 12.27 0 13.22 


To count the fitness for ES algorithm we use Equation (1). 


1 
F=— = (1) 
2, cost, + > penalty. x 1000 
i=l i=l 


2.3. Evolution strategis (ES) 

Evolution Strategy Algorithm using the two main parameters (u and à) where u is the number of 
candidate solutions in the generation of parent while à is the number of candidate solutions produced from 
the parent generation. That’s why it should be better that à value greater than or equal to u (A > p). The 
pseudocode of ES as described by Mahmudy [1] can be seen on Figure 1 and as described by Brownlee [25] 
the flow of ES can be interpreted as Figure 2. 


t=0; 

initialize (P(t=0)); 

evaluate (P(t=0)); 

while isNotTerminated() do 
Pp(t) = selectBest(y, P(t)); 
Pc(t) = reproduce (A, Pp); 
mutate (Pc(t)); 
evaluate(Pc(t)); 
P(t+1) = Pe(t); 
t=ft+l1; 

end 


Figure 1. Pseudocode of ES 
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Initialize Population 


Evaluate Population 


Parent Selection 


Offspring 
Generation 


Survivor 
Selection 


Figure 2. The flow diagram of evolution strategy 


Continue? 


Stop? Termination 


Population will be initialized randomly and will be evaluated using fitness function which will be 
calculated based on the total price of feed ingredients used in feed compositions. Only the best two parent 
will be selected and every parent will be mutated to generate offspring. Finally the offsprings and parents will 
be evaluated, to select which one the best and will be survive. Meanwhile the parameter for ES is use the 
parameter as can be seen on Table 6. This parameter is based on research by [1]. 


Table 6. Parameter for ES 


Generation Number 100 
Mutation Rate 0.6 
Population Size 300 


For this research, there will be a slight modification on the ES algorithm. The modification can be 
seen as red marks on Figure 3. As mentioned by previous study by Sihananto et.al. [1] The purpose of this 
modification is to prevent the existence of genes with negative values. 


t=0; 

initialize (P(t=0)); 

evaluate (P(t=0)); 

while isNotTerminated() do 
Pp(t) = selectBest(u, P(t); 
Pc(t) = reproduce (A, Pp); 
mutate (Pc(t)); 
if gen[i] Of MutatedPc(T)<0 

then gen[i]=random{0,1]; 

evaluate(Pc(t)); 
P(t+1) = Pc(t); 
t=f+l; 

end 


Figure 3. ES modification 
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2.4. Firefly algorithm (FA) 
The pseudocode of Firefly Algorithm as described by Yang [26] can be seen on Figure 4, while flow 
diagram of Firefly Algorithm can be seen on Figure 5. 


Objective function f(x), x= (21,....24)" 
Initialize a population of fireflies x; (i = 1, 2,...,n) 
Define light absorption coefficient 7 
while (t <MaxGeneration) 
for i =1 : n all n fireflies 
for j = 1: i all n fireflies 
Light intensity J; at x; is determined by f(x;) 
if (1; > i) 
Move firefly i towards j in all d dimensions 
end if 
Attractiveness varies with distance r via exp|—r] 
Evaluate new solutions and update light intensity 
end for j 
end for i 
Rank the fireflies and find the current best 
end while 


Postprocess results and visualization 


Figure 4. Firefly algorithm pseudocode [26] 
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Figure 5. The flow diagram of firefly algorithm 
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As usual the population will be initialized randomly, and then among the many fireflies some 
fireflies will be selected as candidate solution. After the cost of the candidate solution firelies have been 
counted, the best candidate will be the ‘center’ and all other fireflies will be ‘fly’ to that center. The penalty 
will be given for each fireflies after each iteration. After maximum iteration, the one with less penalty value 
will be selected as solution. For the parameters, after many trial and errors, we found the optimum parameters 
for Firefly Algorithm (FA) as can be seen on Table 7. To make the test fair, we use 100 iteration both for ES 
and FA. And for every case we repeat the test 10 times because both ES and FA are metaheuristic algorithm 
that may result different solution from time to time. 


Table 7. Parameter for FA 
Iteration 100 
Number of Fireflies 35 


3. RESULTS AND ANALYSIS 
All the testing scenarios are implemented using Java programming language and consist of two 
section: the average fitness test and the running time test. 


3.1. Average fitness test 

The performance comparison of ES and FA can be seen on Table 7. Ther performance is being 
measured by fitness formula as described on (1). The test is being run 10 times to get average fitness as 
shown in Table 8. 


Table 8. Result of Test of Average Fitness 
Average Fitness 


Type of Chicken 


ES FA Pearson's Square 
Layer Starter 0,3306 0,4234 0.01818 
Layer 0,3922 0,4367 0.008 
Layer Finisher 0,3556 0,4369 0.0092 
Broiler Starter 0,3561 0,4406 0.00675 
Broiler Finisher 0,3558 0,4406 0.00675 


As we can see on Table 8., on layer hens case, the chicken growth’s phase that being tested is on 
starter phase, layer phase, and layer finisher phase [23]. The Pearson’s Square got lowest average 
performance among the three methods, while ES got much higher score than Pearson’s Square and FA is the 
highest. 

The pattern on layer hens case repeated on broiler chicken case also can be seen on Table 8., which 
tested on two growth phase: starter and finisher phase [22]. The Pearson’s Square once again got lowest 
average performance, ES got medium score but FA is the highest. 

The ES is perform much better than classic Pearson’s Square method because ES has ability as 
metaheuristic algorithm where it will generate sets of random solutions then evaluate then do selection on 
best solution. This process is called ‘evolution’ by Mezura-Montes et al. [15] and by ‘evolution’ process, the 
better solution is being found. 

But if we compared ES to FA, FA performed better on every test case. This may be happen because 
FA is population based algorithm and social behaviour of its fireflies can be easily formulated and associated 
with given optimization problem [21] and ES like GA sometimes fall into premature convergence, while 
swarm intelligence algorithm like FA can reach wider search area thus prevent the risk of premature 
convergence [27]. 


3.2. Running time test 

The Table 9 show us the comparison of ES and FA performance based on running time in 
microseconds (ms). The PS will not be measured because the PS is being compute by Microsoft Excel and it 
gives almost instaneous result as soon as we change the variable. As shown on Table 9, ES needs more 
execution time than FA. This can be explaining by the characteristics of ES and FA itself. ES has 
characteristic to exploit solutions, it is more a ‘fine tuning’ characteristics than another evolutionary 
algorithm [28], as consequences, it will be slower for ES to reach convergence point. Meanwhile FA is 
powerful on globally-searching promising area. As part of swarm intelligence algorithm, the fireflies of FA 
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can search solution on wider area than ES thus it can find solution faster than ES. Based on this behaviour, 
from the start FA focused on less computational time than ES which need longer computational time to 
achieve result as best as FA. 


Table 9. Result of Running Time Test 


Average Execution Time (ms) 


Type of Chicken ES FA 

Layer Starter 851,174 381,898 
Layer 851,897 374,957 
Layer Finisher 843,925 378,288 
Broiler Starter 819,184 380,966 
Broiler Finisher 831,813 362,346 


4. CONCLUSION 

The implementation of Firefly Algorithm to make chicken feed composition both on layer and 
broiler chicken already proofed that Firefly Algorithm can be used to make chicken feed composition based 
on its average fitness result and computational time. With some forward research, the Firefly Algorithm may 
be implemented as helper algorithm for local farmers on Indonesia, so they can compose their own 
chicken feed. 

Meanwhile Evolution Strategies performance itself is not too bad compared with Firefly Algorithm. 
With some tuning and another modification such as hybridization to improve its computational time and 
fitness performance, the Evolution Strategies itself may become an alternative algorithm to make livestock 
feed composition. 
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